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DESCRIPTION 



Laplace's Equation - Uebmann Method (DOI Only) 

Open 

5 through 46. Grid points and key. words may fill memory 

-8 

Arbitrarily set at 10-5, * u t nay be changed by user. 
3PQ milliseconds, where P » number of interior points 

q » number of times over grid 
For rectangular regions and 5 - pl&ce accuracy Q is 
roughly- approximated by 



where 



K* 

H 

M 



logio K * 

[1 - «*A ([H-l]- 2 + [M-l]" 2 )1 2 , 
number of columns, ° 
number of rows . 



This code gives the solution V(x,y) , corresponding to 
given boundary values, of the difference equation V^ + 
Vw = , which approximates Laplace/s equation 

yy * 

d 2 v + 3 2 v - o 

5x2 Sy 2 

The solution is found for all interior grid points of an 
arbitrary closed region. The scheme used is that of 
Liebmann, which requires an initial approximation for each 
interior point and then improves the approximation by taking 
for a new value at a point the average of the previous 
function values at the four neighboring points of the grid. 
When successive approximations at every point of the grid 
fail to change by more than a predetermined amount, the 
present values of V are taken as the solution. For an 
initial approximation the interior points are taken as 
linear interpolations between theiboundary points marking 
each horizontal sequence of interior points . 
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The Maximum number of boundary and interior points depends on the 
complexity of the boundary but is limited to about 950 far very simple 
boundaries. First, the boundary should be enclosed in a rectangular grid 
of M rows and N columns, preferably with If > M. Known boundary points, 
of course, should correspond to points of the grid. If a given interior 
point has a sequence position i in the MxH grid (reading from left to 
right and top to bottom) then its neighbors have the positions i - N, 
i - 1, i + 1, 1 + N. Since the exterior points will not be stored it is 
necessary to find the actual addresses of the points in sequential positions • 
i, i - H, and i + N. This is done as follows x in preparing the problem 
a count is made of the number of points in each sequence of interior, 
boundary, and exterior points, with the sequences arranged in the same order 
as they are on the grid. These are called temporary key words and are placed 
on the tape after 00 65K as 20 nF 00 F, 30 nF 00 F, ho nF 00 F for 
interior, boundary, exterior sequences, respectively, where n is the number of 
points in a sequence. A tally should be kept of the total number of boundary 
and interior points, for this number plus the number of sequences of interior 
points may not exceed 950. The tape of temporary key words must be terminated 
with 00 F 00 F. It is placed in the memory in locations (65 ff ) later 
occupied by the function values at the grid points. Temporary key words are 
operated on to form the addresses of locations i, i - H, i + H, for each i 
which starts a sequence of interior points, and to place them into one word. 
These words are called the key words and are put into the memory in sequence 
starting at any location following the space required for the 'grid points. 
A suggested position is at address j « [998 minus number of sequences of 
inter ior t points J . Key words are of the form [mx2 ;7 + px2 + q x 2 Jy \ $ 
where N(m), N(p), N(q) are the (i - N)th, i - th, (i + N) - th numbers of the 
grid, respectively. 

Without loss of generality all boundary point values are required 
to lie in the range < V < l/k. Then the interior points are also 
within this range and the boundary point values may be marked by enter leg 
them as negative values . Requiring a maximum of l/k permits adding the 
four points to be averaged without chance of overflow. The boundary ur ; 1 
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values are written as S followed "by up to 11 decimal digits and are placed 
on the tape in sequential order except that after the first boundary point 
following a sequence of n interior points the number X (n + l), where (n + l) 
must be a 3 - decimal digit number, is inserted. This list is terminated 
by the character N. Input of this tape places in order the boundary point 
values aT K* fills in **-li sequences of interior point values by linear inter* 
polations between the boundary values marking each sequence. Boundary 
point values are stored as their negatives* interior points are unchanged. 
Parameters required are (l) N, the number of columns, to be 

placed in location 3 in the form N x 2" 19 , and (2) j, the address for the 

-39 
start of the key words, to be placed in location k in the form j x 2 

.39 
The print parameters are (l) N (3) - k x 2 , where k is the number 

* -.39 

of decimal places to be printed, and (2) H (k) *= s x 2 Jy 9 where s 

is the address of the last point to be printed. 

The complete routine for solving the problem requires the following 

tapes: 

1. Decimal Order Input ( XI - 18 ) 

2. Temporary key words (in 65 ff) and parameters 

3. Routine to form key words (in 10 - 6k), followed by 26 10H 
k. Decimal number sequence input ( N3 - 23 ) and modification 

(in 26 - k6). N3 is obsolete, but it has been copied into 
Gl ftT ^ no copy is required by the programmer. 

5. Boundary input routine (in 5 - 25), followed by 26 28JT 

6. Boundary points, followed by N 

7. Parameters again and routine evaluating interior points 
(in 10 - 38), followed by 2k ION 

8. Print parameters, print routine, and Single Column print 
( P6 - 25 ) (in 10 - 28), followed by 2k ION 

The order listed must be adhered to although only 2, 6, and parameters need 
to be constructed for each problem. Codes k and 5 overwrite 3* 7 overwrites 
k and 5, and 8 overwrites 7« 

Relative location 28 in the code for evaluating interior points 
(fixed location 38) contains the constant which determines the accuracy 
and thus the stopping point of the routine . This has been Bet arbitrarily 



Jt. 



at 5 x 10 to give five-place accuracy, but say be varied at the 
discretion of the user by setting I(38f) ■ 5 x 10"( a + *), (vnere m la 
the number of declaal places of accuracy desired) after input of the routine 
for evaluating interior points, hut before its execution. 



An example follows: 



M*5 



ti^io 




Suppose the encircled boundary points have the value 10 and the others the 
value (-2). First add 2 to the function to get it non-negative through- 
out the region. Then scale by a factor of 50 so that the function values 
V are all in the range ( < V <,2k ) . The temporary keyB are kO IF 00 T, 
30 7F 00 F, kO 2F 00 F, 50 IF 00 F, 20 7F 00 F, 50 IF 00 T, kO IF 00 T, 

50 IF 00 F, 20 5F 00 F, 50 5F 00 F, 20 2F 00 F, 50 2F 00 F, 20 2F 00 F, 

50 IF 00 T, 1*0 3F 00 F, 50 IF 00 F, 20 IF 00 F, 50 U 00 F, Ko IF 00 F, 

50 2F 00 F, kO 5F 00 F, 50 IF 00 F, 00 F 00 F. Parameters are U(5) « 10 x 2" 19 

and N(4) « [998 - 5] x 2~ 39 = 995 x 2~ 59 . The list of boundary points is 32k, 

S24, 32k, S, S, S, S, S2k, 3, K008, 32k, 32k, KX*, S, S, S, K003, 32k, 32k, 

K005, S, S, K002, 32k, 32k, 3, 5. 



Igr 



October 



MSB 

CQDBD BYj 
AiraOVKD BT 



SL 



&J»» ■PlQ/*. 



lf«M.5 



;CIOW 



ORDER 



NOTES 



1 G 1 



3 
4 

5 
6 

7 



Q 



10 



II 



13 
j.ij. 

'-5 
16 



00 10K 

41 7P 

40 8F 

50 (65 )F 

01 3F 
LO 52L 
36 4L 

41 ( )F 
24 999F 
LO 52L 
LO 52L 
36 47L 
41 5F 
40 6f 

L5 7F 
LO 8F 
26 8l 
40 F 
L5 1L 
46 26L 
46 42L 
LO 51L 
h6 11L 
50 ( )F 
01 3F 
LO 53L 
36 18L 
11 3F 
L5 5F 
S4 F 
40 5F 
LO 3F 
32 20L 
L5 11L 
LO 51L 



1(7) - H(8) « 



From 46 r j— Q « current temporary key 



No 
By 41 

From 2* 

No 
Yes 



_ Test for end of temporary keys 

Insert at end of keys 
Read in "boundary input 11 
Test for interior sequence » 



1 

l i 

I L 



By 10* 17 p 
From 17* 

Yes 

No 
From 20 



Yes 
No 



N(5) «= N(6) « 

N(0) » [number of interior parts 

4- number of boundary points] 2""" 1 "* 



Plant address- of current temporary key 



plant address of previous teasporsury keif 



Test for exterior a aquas 



N(5) = [number of points we nave racked} 
up] • 2" 19 I 

Have we backed up a row? 



Plant address of previous 



-ry fee^l 



■ --Tcgr 


ORDER 


■0TH3 PAGE 2 CI 




46 UL 








26 11L 






lb 


11 5F 


From 12» 






L5 6F 
S4 F 




H(6) « [number of exterior point® 

-19 
we hare backed up] . 2 




40 6? 




_ 


20 


22 13L 
L5 54L 


From 15* 


-19 
A ■ [address of first mesh word] . 2 


21 


L4 f 






22 


40 IF 
LO 3F 
L4 6f 




N(l) * m . 2" 9 + p . 2~ 19 


23 


00 10F 
L4 IF 








40 IF 




« 




41 5F. 






?5 


4l 6F 




N(5) * N(6) » 




26 26L 




. 


26 


50 ( )F 


By 9,p 






01 5F 


From 32' 


Test for exterior sequence 


27 


LO 53L 


Yes 






56 33L 






;.3 


11 3F 
L5-5F 


No 
From 35 


" N(5) = [number of points we have 

-19 
progressed . 2 ^ 


jj 


S4 F 




• 


•:; -i 


40 5F 
LO 3F 




Have we progressed a row? ,' 4 , 




32 35L 


Yes 


aJ 


M 


L5 26L 


No 






. L4 51L 




Plant address of next temporary key 


52 


46 26L 








26 26L 


• " 


■* 


55 


11 3F 


From 27' 




5 1 * 


L5 6F 
S4 F' 




N(6) = [number of exterior parts we 

-19 t 
have progressed] . 2 




40 6F 




I 



LOCATION 



i.-~ 



,35 
56 

37 
.58 

39 

14-0 

in 

42 

43 
44 

45 

46 



48 
52 



ORDER 



22 28L 
L5 IF 
L4 3F 
LO 6F 

10 20F 
42 IF 
L5 IF 
26 39L 
40 ( )S4 
L5 39L 
L4 51L 
46 39L 
46 3L 
26 42L 

50 ( If 

00 3F 

11 3F 
L5 7F 
S4 F 
40 7F 
L5 1L 
L4 51L 
46 1L 
26 1L 
LO 52L 
32 48L 
26 43L 
11 3F 
L5 8F 
S4 F 
40 8F 
22 43L 
00 IF 
00 F 
00 F 
00 IF 



HOTES 
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G 1 I 



From 30 



By 9* 

From 48 
From 50' 



From 5 
No 
Yes 
From 47' 



u 



,-9 



H(l) = m . 2 y + p . 2" x ' f q • 2 ' 



19 



~*Q 



__ = key word 



_ Store key word 



_ Plant addresses for next key word 



Q = [number of parts in present interior! 
sequence J . 2 j 



N(7) = [number of interior parts + 



number of exterior points + number of * 



boundary points] . 2 



-19 



Set address for next temporary key 



Test for boundary sequence 



N(8) = [number of exterior points] x 



2 -19 

0-39 



*ETTjr5wa»jflusfM6iieaai 



ORDER 



00 F 
00 4F 
jk 00 47F 

00 F 

26 ION 
BOUNDARY INPUT 



56 



37 



58 



39 



1*0 



00 26K 

DNSI Code N3 

00 36K 

66 IF 

2.6 5F 

40 (kj)F 

L5 37F 

L4 32F 

46 57F 

L5 2F 

LO 45F 

40 IF 

54 999F 

00 5K 

55 IF 
52 1L 
26 57F 
40 IF 
L5 57F 
LO 20L 
k6 TL 
46 10E 
46 l4L 
LO 20L 
46 6L 
46 13L 
L5 ( )F 
50 21L 

LO ( )F 
40 4F 



NOTES 
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By 5 
By 3 



4 x 2"^ 
4? x 2~ 19 



26 - 46F 



Modification to DNSI code (N3) 



Read in "Evaluate interior points" routine i 



interior points or boundary points 

N(i) = 10~ 5 (n + 1) 

Get address of inter icr points 



Plant first interior point address 



__ Plant left hand boundary point address 



Find functional increment per interior 
point 



| LOCATION 


ORDER 




NOTES PAGE 5 


G 1 


8 


75 ^F 
66 IF 










9 


S5 F 




■ 


» ■ 




10 


L5 ( )F 
kO 3F 


By 3' 


■ 


N(3) = right hand boundary point 




11 


L5 21L 
L2 IF * 






1 Count to n and test 




12 


^0 II. 
36 18L- 










13 


L7 ( )F 

Ik hF 


By 5 ! 


,15" 


Form initial values for interior 


points 


Ik 


hO ( )F 
L5 l^L 


By 4, 


16. 


■ 




15 


\6 15L 
lA 20L. 










16 


h6 ikL 
00 63F- 






Set addresses in loop 




11 


00 63F 
26 11L. 










.18 


L5 ikL, 








1 




h6 3TL 






Prepare to store right hand bound 


stry 


19 


L5 3F" 
26 37F 






point 




20 


00 IF 
/I 00 1F ; 






2 =19 + 2 -59 




21 


1 « DOFJ 00001 
1 1 00000 1000J 
// 26 28N 






lo~ 5 + io° 9 














EVAI 


,UATE INTERIOR PO 


INTS- 










00 10K" 













hi 6f 






_N(6) - 






L3 ( )Sk 


By 23 




Test for end of grid 




1 


32 19L 




. 


1 






L5 ( )S* 


By 23' 


_____ 


N(2) = current key word 





WCMEXM 


UKIUBt 


■DSB K^ 


2 


fcO 2P 


J 


1 




k6TL 






3 


MlkL 




_ Plant p, q. 


! 


k2 11L 






k 


IA 25L 


■> 






1*6 13L 


_J° Plant p «- 1 


5 


LO 26L 
1*6 12L 


T 

J Plant p - 1 


* 


10 50F 








U2 12L 




Plant m 


7 


L5 (p)F 
1*0 4f 


By 2« 1 M(p) « 1(4) 


8 


36 111 




Test far boundary point 




L5 L 


Yes 


♦ 


i ^ 


lA 25L 


• 


Hant key work addresses 


| 


llS L 




» 


1 10 


1*2 1L 








22 L 




Sad of row 


i 

f U 


19 37F 


From 8 




j 


L7 (Of 


3y 3' 




j 12 


L6 (p-l)F 


By 5 f 


H i+1 (p) « lA tl i+1 (a) * It(q) 


1 


US (a)F 


3y 6' 


_ (p-1) + I^P+l)] 


13 


L6 (p+l)F 
10 2F 


By V 




Ik 


VO (p)F 








LO k¥ 






15 


kO kT 








L7 ^F 




.1(6) -E |h 1+1 (p) - S t (p)i 


16 


lA 6F 








1*0 6f 






17 


36 18L 




Test for overflow 




1*0 27L 


Yes 


H(2fL) < 


18 


L5 2F 






» 


lA 25L 




Step he, p, q. by 1 


^ 


26 2L 








12 27L 


Frea 1 


mm 



6 



mm 



JUOCATION 



ORDER 



20 56 22L 
L5 28L 

21 L0 6F 
ik 999F 

22 ^9 27L 
L5 2^L 

25 ^2 L 
k2 1L 

24 26 L 
OF SV 

25 00 1025F 
00 1025F 

26 00 1026F ' 
00 1025F 

27 40 F . 
00 F 

28 f °° F 0000 ° 
005 000000J 

2k ION 

PRINT ROUTINE 

i, 00 10K 

L5 (V7)F 
50 L 
26 16F 
19 18F 
I> L 
46'L 
LO 5L 
52 ^L 
26 L 
00 F 
L5 1SU 
50 L 
00 16K 



Ho 



From 20 



U 



HOSES 



Test for overflow 
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Test for end 
J Read in print routines 
3- Reset H(27L) > 






Reset addresses and start over 01 

Address used 

2-9 + 2- 1 9 + 2- 2 9 + 2-39 



2~ 9 + 2 . 2-^ + 2- 2 9 



1/2 



5 . 10' 






Argument to "single colims 



Step address for new 

M Test for end 

Print ngxx point 
final stop 
End constant 



Single Column Print 
2k ION 



(P6) 

_L 



